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ABSTRACT 

We discuss the evolution of the magnetic field of an accreting white dwarf. We calculate 
the ohmic decay modes for accreting white dwarfs, whose interiors are maintained 
in a liquid state by compressional heating. We show that the lowest order ohmic 
decay time is (8-12) billion years for a dipole field, and (4-6) billion years for a 
quadrupole field. We then compare the timescales for ohmic diffusion and accretion at 
different depths in the star, and for a simplified field structure and assuming spherical 
accretion, study the time-dependent evolution of the global magnetic field at different 
accretion rates. We neglect mass loss by classical nova explosions and assume the white 
dwarf mass increases with time. In this case, the field structure in the outer layers of 
the white dwarf is significantly modified for accretion rates above the critical rate 
M c ~ (1-5) x 10~ 10 Mq yr _1 . We consider the implications of our results for observed 
systems. We propose that accretion-induced magnetic field changes are the missing 
evolutionary link between AM Her systems and intermediate polars. The shorter ohmic 
decay time for accreting white dwarfs provides a partial explanation of the lack of 
accreting systems with r* 10 9 G fields. In rapidly accreting systems such as supersoft X- 
ray sources, amplification of internal fields by compression may be important for Type 
la supernova ignition and explosion. Finally, spreading matter in the polar cap may 
induce complexity in the surface magnetic field, and explain why the more strongly 
accreting pole in AM Her systems has a weaker field. We conclude with speculations 
about the field evolution when classical nova explosions cause the white dwarf mass 
to decrease with time. 
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1 INTRODUCTION 

There are now about 65 isolated magnetic white dwarfs 
known, and a further 90 accreting magnetic white dwarfs 
(see Wickramasinghe & Ferrario 2000, hereafter WF, for 
a recent review). Interesting differences are emerging be- 
tween the isolated and accreting populations. Measured field 
strengths of the isolated WDs range from 3 x 10* G to 10 9 G, 
whereas the range of field strengths of accreting WDs may be 
smaller. The magnetic field strengths in the AM Her binaries 
are directly measured to be in the range 10 7 -2 x 10 8 G, and 
the magnetic fields of the intermediate polars (IPs) are in- 
ferred to range from 10 7 G down to « 10 5 G. The observed 
fraction of magnetic systems is 5% for isolated WDs, but 
25% for accreting systems. There is evidence that isolated 
magnetic white dwarfs are more massive than the typical 
0.6 Mq non-magnetic WD, with a mean mass >0.95 Mq. 
Measurements of accreting magnetic white dwarfs, however, 
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give masses ~ 0.7 Mq, consistent with measured masses of 
non-magnetic (B < 10 5 G) accretors. One similarity between 
isolated and accreting WDs is that the magnetic field is of- 
ten seen to be complex, with a significant quadrupole (or 
offset dipole) component. 

The origin and evolution of white dwarf magnetism is 
an important, but as yet unsolved, part of our understand- 
ing of stellar evolution. Several papers have discussed the 
evolution of magnetic fields in isolated, cooling white dwarfs. 
Ohmic decay of magnetic fields was first addressed by Chan- 
mugam & Gabriel (1972), and Fontaine, Thomas, & van 
Horn (1973), who showed that the decay time of the low- 
est order decay eigenmodes was longer than the evolution- 
ary time of the white dwarf. Wendell, van Horn, & Sar- 
gent (1987; hereafter WVS) followed the evolution of dif- 
ferent decay modes along a cooling sequence. They found 
that the fundamental dipole mode decayed by a factor of 
two in 10 10 years. These studies suggested that WD mag- 
netic fields are fossil fields left over from previous evolution. 
For example, one proposal consistent with the high masses 
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(>0.95 Mq) of magnetic white dwarfs is that they result 
from evolution of the strongly magnetic (10 2 -10 4 G) Ap and 
Bp main sequence stars. The higher order field components 
decay more rapidly, however, leading Muslimov, van Horn, 
& Wood (1995, hereafter MVW) to study the Hall effect as 
a way to generate field complexity during the lifetime of the 
white dwarf. 

The result that white dwarf magnetic fields evolve only 
very slightly over their lifetime has been applied in almost all 
studies of accreting white dwarfs. For example, in evolution- 
ary studies of magnetic cataclysmic variables, it is presumed 
that the magnetic field of the white dwarf does not change 
with time. One exception is King (1985), who applied the 
results of Moss (1979) to accreting systems, and suggested 
that meridional currents could submerge magnetic flux in 
the outer layers of rotating magnetic white dwarfs. It is the 
purpose of this paper to show that the process of accretion 
itself may significantly change the surface magnetic field of 
an accreting white dwarf. 

Many authors have studied the effects of accretion on 
the magnetic field of accreting neutron stars. This is moti- 
vated by the idea that the rapidly rotating millisecond ra- 
dio pulsars are produced by accretion, which spins up the 
neutron star to short periods, and perhaps causes a reduc- 
tion in magnetic field strength from the 10 12 G fields seen 
in most radio pulsars to the 10 8 -10 9 G fields of the mil- 
lisecond pulsars (see Bhattacharya 1995 for a review). In a 
recent paper, Cumming, Zweibel, & Bildsten (2001; here- 
after CZB) returned to the suggestion that the accretion 
flow directly screens the internal magnetic field (Romani 
1990, 1995). They compared the timescales for ohmic dif- 
fusion and accretion in the thin outer layers of the neutron 
star, and computed steady-state magnetic profiles, taking 
the field to be horizontal and the accreted matter to be un- 
magnetized. For this simplified geometry, the field was found 
to be strongly screened for accretion rates greater than the 
critical rate M ~ 10~ 10 Mq yr" 1 , whereas at lower accre- 
tion rates, the underlying magnetic field was able to pen- 
etrate the freshly accreted material. This result fits nicely 
with the observation that the only weakly-magnetic accret- 
ing neutron star to show its magnetic field directly (the ac- 
creting millisecond X-ray pulsar, SAX J1808. 4-3658; Wij- 
nands & van der Klis 1998; Chakrabarty & Morgan 1998) 
has a lower accretion rate than other systems (time-averaged 
M w 10 -11 Mq yr" 1 ), low enough that screening would not 
be effective in hiding its field. 

A crucial question for the evolution of accreting white 
dwarfs is whether the white dwarf mass increases or de- 
creases with time. For accretion rates < 10 -7 Mq yr" 1 , the 
accreted hydrogen and helium burns unstably, giving rise to 
classical nova explosions (e.g., see Fujimoto 1982; MacDon- 
ald 1983). The observed heavy element enrichment of nova 
ejecta has been used to argue that the nova explosion ex- 
cavates material from the white dwarf, so that its mass is 
decreasing with time (e.g., Livio & Truran 1992). However, 
this is a somewhat open question, both because of uncer- 
tainty in measurements of ejecta masses, and in theoretical 
modelling. Theoretical calculations find that mixing of extra 
CNO nuclei into the burning shell is needed to achieve rapid 
enough energy release to drive novae with rapid rise times. 
However, the mixing mechanism is presently unknown. Pri- 
alnik and Kovetz (1995) addressed this question with hy- 



drodynamic simulations including diffusion, and found that 
the white dwarf mass decreased for M < 10~ 9 Mq yr _1 , 
and increased for M > 10~ 7 Mq yr -1 . At intermediate ac- 
cretion rates, the mass loss depended on the white dwarf 
temperature. The accretion rates of IPs are estimated to be 
~ 10 -9 Mq yr~ x (see discussion in §5.1), and their core tem- 
peratures in the middle of the range considered by Prialnik 
and Kovetz (1995) (see §3.2). On the other hand, AM Her 
systems accrete at typical rates « 5 x 10~ n Mq yr -1 (§5.1). 
Thus a range of behaviour might be expected within the 
population of accreting magnetic white dwarfs. The mag- 
netic field of the white dwarf itself may affect the nova out- 
burst, for example by enhancing mass loss in the equatorial 
plane (Livio, Shankar, & Truran 1988; Livio 1995). 

In this paper, we assume that the white dwarf mass 
increases with time during accretion, and study the conse- 
quences for the evolution of the white dwarf magnetic field. 
We leave the case of decreasing white dwarf mass for a future 
study. We first show that comparing accretion and ohmic dif- 
fusion times in white dwarfs gives the same critical accretion 
rate M c w (1-5) x 10" 10 M yr 1 as found for neutron stars 
by CZB (for reasons we describe in §3), and then go on to 
calculate the evolution of the magnetic field. We show that 
the resulting field evolution may explain several properties 
of observed systems. 

A complication in calculating the evolution of the mag- 
netic field in both accreting neutron stars and white dwarfs 
is that the accretion flow is channelled onto the magnetic po- 
lar cap (for white dwarfs, this occurs for B > 10 G). Rather 
than tackle the complex problem of the subsequent spread- 
ing of matter and evolution of the magnetic field, CZB pre- 
sumed that the neutron star magnetic field was buried by 
the accretion flow into a flattened configuration, and then 
used a plane-parallel model to ask what accretion rate was 
needed to keep the field from reemerging by ohmic diffusion. 
In this paper, the simplification we make is to study the evo- 
lution of the global magnetic field under spherical accretion. 
The accreted matter is expected to spread away from the po- 
lar cap when the magnetic tension force no longer supports 
the hydrostatic pressure (e.g., Hameury et al. 1983), so that 
spherical accretion is a good approximation except in a thin 
layer of mass AM < 10" 10 Mq B 2 for a 0.6M Q white dwarf, 
where B7 — B/W 7 G is the surface field strength. We in- 
clude the uncertainty of what happens in this thin spreading 
layer as an uncertainty in the surface boundary condition for 
the spherically accreting models (see discussion in §4.1). 

We start in §2 by considering ohmic decay in accreting 
white dwarfs, which have liquid interiors because of compres- 
sional heating by accretion. In §3, we compare the timescales 
for ohmic decay and accretion, and relate our results to those 
of CZB. In §4, we calculate the global evolution of the field, 
under simplifying assumptions about its geometry, for dif- 
ferent accretion rates. We relate our models to the recent 
calculations of Choudhuri & Konar (2002) for accreting neu- 
tron stars, which adopt a similar approach. We find that for 
M > M c , the surface field strength is reduced by accre- 
tion because ohmic diffusion is not rapid enough to allow 
penetration of field into the newly accreted layers. In §5 we 
discuss the implications of our results for observed systems. 
In §6, we conclude with some speculations about magnetic 
field evolution in the alternative scenario of decreasing white 
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dwarf mass with time. In the Appendix, we describe our cal- 
culations of electrical conductivity for arbitrary degeneracy. 



2 OHMIC DECAY IN LIQUID WHITE 
DWARFS 

In this section, we discuss ohmic decay of the magnetic field, 
neglecting the effects of accretion for now. In an accreting 
white dwarf, compressional heating of the core results in a 
central temperature > 10 7 K (Nomoto 1982), so that the 
interior is kept in a liquid state. In that case, the electrical 
conductivity is independent of temperature, allowing us to 
use simple zero-temperature white dwarf models. We calcu- 
late the ohmic decay modes following the original work of 
Chanmugam & Gabriel (1972), Fontaine et al. (1973), and 
WVS. We show that the lowest order decay mode has a de- 
cay time w (8 — 12) x 10 9 yrs, depending on white dwarf 
mass. 



2.1 Electrical conductivity and decay time 

Ohmic decay is described by the induction equation 

f = -Vx ivVxB), (1) 

where r\ = c 2 /4tvu is the magnetic diffusivity, and a the 
electrical conductivity. The associated characteristic decay 
timescale is « 4iYaL 2 /c 2 , where L the lengthscale over which 
B changes. 

Compressional heating maintains the interior of most 
accreting white dwarfs in a liquid state. The melting tem- 
perature T me it is determined by the condition F « 173 (e.g., 
Farouki & Hamaguchi 1993), where F = (Ze) 2 /ksTa mea- 
sures the ratio of Coulomb to thermal energy. For simplicity, 
we write F for a single species of ion with charge Z, mass A 
and interion spacing a, giving 



T melt ~ 3 x 10 6 K ^ (f ) ) • (2) 

Compressional heating gives core temperatures > 10 7 K > 
Tmcit (Nomoto 1982), thus we expect most accreting white 
dwarfs to have liquid interiors. It is possible that some mas- 
sive white dwarfs with central densities > 10 s g cm -3 ac- 
creting at low enough rates have solid cores, but we do not 
consider these systems here. 

The electrical conductivity in the liquid interior is set 
by collisions between the degenerate electrons and the non- 
degenerate ions. The conductivity is 



8.5 x 1(T 



A ei {Z) 



1+x 2 



(3) 



(e.g. Yakovlev & Urpin 1980; Itoh et al. 1983; Schatz et 
al. 1999), where x = pF/m e c measures the Fermi momentum 
Pf of the electrons, A e i w 1 is the Coulomb logarithm, and 
(Z) is a measure of the mean nuclear charge, given by (Z) = 
fi e ^2 XiZf /At, with Zi the nuclear charge, Ai the mass, and 
Xi the mass fraction of species i, and p, e the mean molecular 
weight per electron. In terms of x, the density is p = 1.9 x 
10 6 g cm" 3 x 3 (Me/2). 

The fact that conductivity is independent of tempera- 
ture allows us to adopt simple zero-temperature white dwarf 



models. In the outer envelope, where the electrons are non- 
degenerate, the conductivity is temperature-dependent. In 
the Appendix, we describe a method for calculating the elec- 
trical conductivity for arbitrary degeneracy. However, the 
results of this section are insensitive to the conductivity in 
the outermost layers. We adopt pressure as the independent 
variable, and integrate the equations of hydrostatic balance 
and mass conservation outwards from the centre of the white 
dwarf for an initial choice of central density. We stop the in- 
tegration once the pressure falls below some fixed fraction 
of the central pressure, at which point the radius and mass 
are well-determined. We assume the composition is an equal 
mixture by mass of carbon and oxygen, giving fi e — 2 and 
(Z) = 7. 

Some properties of these zero-temperature white dwarfs 
are shown in Figure [lj. We show the radius, central den- 
sity, central conductivity cr c , and ohmic time at the centre, 
£ohm = AtkR 2 <j c j c 2 , as a function of white dwarf mass. 

Despite the large variations in R and a c , the ohmic 
time at the centre does not depend very strongly on white 
dwarf mass, ranging from (2-6) x 10 11 yrs. To see why, we 
note that at the centre of the white dwarf x> 1, so that 
a oc x oc /5 1//3 . However, if the mass-radius relation is R oc 
M" 1/3 , then p oc M 2 , giving a oc M 2/3 . The ohmic time 
tohm oc R 2 a is thus approximately constant with mass. As 
the white dwarf mass increases, the increase in conductivity 
is offset by the decreasing radius. Most of the variation in 
tohm seen in Figure [l] comes from the Coulomb logarithm 
(see Appendix for our calculation of A e i). 

The difference between ohmic decay times in isolated 
white dwarfs with liquid and solid cores has been pointed 
out previously, particularly by MVW. For a solid, the con- 
ductivity is set by collisions of electrons with phonons, which 
become more common with increasing temperature (for ex- 
ample, see Yakovlev & Urpin 1980; Baiko & Yakovlev 1995). 
An example is Figure 2 of WVS, in which both the conduc- 
tivity and decay timescale are roughly constant for the first 
~ 3 x 10 s years of cooling history when the core is liquid, 
but then increase with time as the core becomes solid and 
cools. 



2.2 Ohmic decay modes 

The ohmic time at the centre t hm is an overestimate of the 
time for magnetic field decay, since both the conductivity 
and radial lengthscale are largest at the centre. To do better, 
we now find the ohmic decay eigenmodes for liquid white 
dwarfs, following the original work of Chanmugam & Gabriel 
(1972), Fontaine et al. (1973), and WVS. 

We assume an axisymmetric poloidal magnetic field, 
B(r,9,t) = B r (r,8,t)e r + Be(r,9,t)ee which we write in 
terms of a vector potential B = Vx A where A — 
A^ir, 9,t)e$ (see Mestel 1999 for a useful discussion). Us- 
ing the identity Vx {e^/rsmff) — 0, we find 



■xV(rA^ sine) , 



(4) 



so that the quantity rA^ sin 9 labels the magnetic field lines. 

Rewriting equation (|jj) , we find the time evolution of A 
is given by 

dA 

— = -i)VxVxA. (5) 
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Figure 1. Properties of the zero-temperature white dwarf mod- 
els. We show the radius, central density, conductivity, and ohmic 
time t h m = AirR 2 a c /c 2 as a function of mass. The conductivity 
is calculated assuming a liquid interior. 



We now seperate into radial and angular pieces, 



(6) 



where P; 1 is the associated Legendre function of order 1. The 
magnetic field components are 

1(1 + 1) 



and 



= -/ ^ i p i (cosr?), 



(7) 



(8) 



where Pi is the Legendre polynomial. We use spherically 
symmetric white dwarf modelsn in which the diffusivity de- 
pends on radius r only, and not on latitude 9. In this case, 
substituting equation (g) into equation (^) gives 

1 Except in the very outermost layers, it is always a good ap- 
proximation that the hydrostatic pressure P is much larger than 
the magnetic pressure B 2 /8n. 
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Figure 2. Ohmic decay modes for a 0.6Mq white dwarf, (i) In the 
upper panel, we show the conductivity as a function of radius. The 
central value of conductivity in this model is Rj 10 21 s — 1 , dropping 
to close to the Spitzer value fa 5 X 10 17 s _1 (eq. Jl9| ) at the outer 
boundary, (ii) In the middle panel we show the n = 1 and 2 
dipole (I = 1) decay modes. We normalize each eigenfunction so 
that its maximum value is unity. The nth mode has (n — 1) nodes. 
The dashed curve shows the n = 1 mode for constant conductivity 
with radius, (iii) The lower panel shows the distribution of current 
density for the eigenmodes shown in the middle panel. 



dRi 
dt 



l(l + l)Ri 



(9) 



a diffusion equation for the radial function ii; . 

Neglecting mass and density changes due to accretion, 
the conductivity of a liquid white dwarf is independent of 
time, and depends only on position, r)(r). In this case, we 
follow WVS, and look for a solution which decays exponen- 
tially with decay time n n , 



Ri{r,t) = CinX tn (r) exp(-t/ri„) . 



(10) 



where C; n is a constant giving the contribution of mode n, 
and the differential equation 
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Table 1. Ohmic decay times (tj„ in units of 10 9 years) 



M = 0.6 M 
(tohm = 3.4 X 10 11 yrs) 
n = 1 n = 2 n = 3 
2 = 1 7.6 2.5 1.2 
1 = 2 3.5 1.6 0.87 
Z = 3 2.0 1.1 0.66 



M = 1.0 M 
(tohm = 5.4 X 10 11 yrs) 
n=l n = 2 n=3 
12 4.1 2.0 

5.7 2.5 1.4 
3.3 1.7 1.1 



d 2 X ln 1(1 + i)X tr . 



dr 2 



+ 



Xi, 



= 



(11) 



determines the eigenfunction Xi n (r). 

We adopt the following boundary conditions (see dis- 
cussion in §5.6 of Mestel 1999). Inspection of equation ( |i~ij ) 
shows that a finite solution at the centre of the star must 
have Xi n oc r l+1 , giving 

dXm (I + l)X ln 



dr 



0. 



(12) 



Outside the star, where the current J oc V X B = 0, we 
choose the solution which remains finite for large r, Ri n oc 



Thus, we take 

IX l n 



R 



(13) 



dr 

at the surface. 

Solution of equation ([H]) requires the conductivity as a 
function of radius. This is shown in the top panel of Figure 
|i| for a O.6A/0 white dwarf, normalized to the central value, 
a c = 9.6 x 10 20 s _1 . The resulting n = 1 and n = 2 dipole 
(I = 1) decay modes Xi n are shown as solid lines in the 
middle panel of Figure ^. The lower panel shows the current 
density J = (c/47r)Vxi?. Using equations (5), (6), and 
(10), this may be written in terms of Xi n as J — e$, 
where 



■hi, 



o (X ln \ 
cn n v r J 



(14) 



Note that the nth mode has n — 1 nodes, the current revers- 
ing direction at a node. 

Table |l| gives the decay time for modes with I = 1,2 
and 3, and n = 1,2, and 3 for 0.6 Mq and 1.0 Mq white 
dwarfs. For both white dwarf masses, we find that the lowest 
order decay mode (1 = 1 and n = 1) has a decay time 
T11 ~ t o hm/40, roughly a factor of four smaller than the 
analytic result for r\ independent of radius (in which case 
Tu = t ohm /iv 2 , see WVS eq. [10]). 

The decay times in Table |l| agree well with the calcula- 
tions of WVS for liquid white dwarfs (compare their Figure 
2). However, we do not agree with the later calculations 
of MVW, who report decay times an order of magnitude 
smaller, (3 — 10) x 10 8 yrs, for liquid white dwarfs (see their 
§3.3), despite using the same conductive opacities for the 
liquid interior as WVS (see their §2.3). The origin of this 
difference is not clear. The central conductivity quoted by 
MVW for a 0.6 Mq white dwarf (« 2 x 10 21 s" 1 ) agrees to a 
factor of two with the value in both our models and those of 
WVS. In addition, changing the conductivity profile seems 
unlikely to give an order of magnitude change in the decay 
time. The difference may be due to an incorrect normaliza- 
tion used by MVW (H. van Horn, private communication). 

To test the influence of the conductivity profile, we com- 
pared our eigenfunctions with those of WVS and MVW (Fig- 



ure 3, left panel, of WVS; Figure 1, top panel of MVW). The 
agreement is good for times > 10 9 yrs, but differs for earlier 
times, for which the WVS and MVW eigenfunctions peak at 
larger radius (r/R « 0.8 rather than r/R sa 0.6). We have 
been unable to identify the reason for this difference. The 
peak at larger radius likely reflects a shallower conductiv- 
ity profile. For example, in the extreme case of constant 
conductivity with radius, the currents are less centrally- 
concentrated than in realistic models. The reason for this 
is that the currents flow in the regions of high conductivity; 
this region is more extended in the constant conductivity 
model. The dashed lines in the middle and lower panels of 
Figure || show the n = 1, I = 1 eigenmode for constant con- 
ductivity (in this case X11 oc rji(r), where ji is the spherical 
Bessel function of degree 1, see WVS). As noted above, the 
lowest order decay time for constant conductivity is a factor 
»; 4 longer than the realistic white dwarf models. 

In summary, we find that the ohmic decay time lies 
in the range ~ (8-12) billion years for a dipole field, and 
~ (4-6) billion years for a quadrupole field, depending only 
slightly on white dwarf mass. 



3 ACCRETION VS. DIFFUSION TIMES IN 
THE ENVELOPE 

The long timescale for ohmic decay raises the question of the 
effect of the accretion flow on the magnetic field structure. 
The decay timescale of « 10 10 yrs is the time to accrete 
the whole star at 10 -10 Mq yr _1 . For more rapid accretion 
than this, one might expect the magnetic field structure to 
be significantly altered as the magnetic field is advected by 
the accretion flow. 

In fact, the evolution of the magnetic field depends 
on the local accretion and diffusion timescales at different 
depths in the white dwarf envelope. In this section, we com- 
pute these timescales, before moving on to simple models of 
the global field structure in §4. 



3.1 Analytic estimates 

We consider the outer layers of the star, where the grav- 
ity g = GM/R 2 is constant. We define the column depth 
y = — J pdz = AM/4-kR 2 , where AM is the mass above 
column depth y. Hydrostatic balance dP/dz = —pg be- 
comes dP/dy = g, giving P = gy. The time to accrete the 
matter above a given column depth is t accr = y/rh, where 
rh = M/4-kR 2 is the local accretion rate per unit area. We 
compare this with the ohmic diffusion time across a scale 
height, t<na = 4tv(tH 2 /c 2 . 

When the electrons are degenerate, the conduc- 
tivity is given by equation (|^). In the outer layers, 
the electrons are non-relativistic, so that P = 6.8 x 
10 20 erg cm -3 p^ 3 (2/ /-i e ) 5 ^ 3 ■ The pressure scale height is 



H = P/pg = 6.8 x 10 7 cm (p 2/3 /g 8 )(2/ Me ) 5 / 3 , 



giving 



tdift = 8.9 x 10 s yrs 



7/3 



A« (Z) g 2 



13/3 



(15) 



To evaluate the accretion timescale, we write the local accre- 
tion rate rh = rh_3l0 -3 g cm -2 s _1 in terms of the global 
rate M = M_i 10" 10 Mq yr" 1 , as rh = M/4tt_R 2 , giving 
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m-3 



0.51 



M-ii 
Rl 



(16) 



If accretion is not occuring over the whole surface of the 
white dwarf, for example magnetically-controlled accretion 
onto the polar caps, the relevant quantity is the local ac- 
cretion rate; we work in terms of the global rate only for 
convenience. We find 

5/3 

(17) 



facer = 4.3 x 10 8 yrs 



Ha 



gsM- 10 

The ratio of tdie to i accI depends on the quantity gR 2 , 
which we rewrite in terms of the mass of the star, gsRy = 
1.3(M/M Q ). This gives 



fdiff 
^accr 



1.6 



2/3 
P 5 



A ei (Z) 



/M©\ 
\ M ) 



8/3 



(18) 



which depends only weakly on depth. 

Equation ([H^) shows that teas ~ facer for accretion rates 
of order 10" 10 M yr" 1 . This is roughly the same critical 
accretion rate as found for the degenerate ocean of an accret- 
ing neutron star in CZB (compare eq. [29] of that paper). 
The reason is that fdiff/faccr oc 1/gR 2 , so that it depends 
only on the mass M and not separately on the gravity and 
radius. 

At densities < 10 3 g cm -3 , the electrons become non- 
degenerate, in which case the conductivity is given by 
Spitzer's formula (Spitzer 1962) 



7.5 x 10 1 



3/2 



A ei (Z) ' 



(19) 



where the temperature T = T7IO 7 K. For an ideal gas, the 
pressure scale height is H = P/pg = kpT '/ pm p g, where p, is 
the mean molecular weight, giving H = 8.3 x 10 6 cm Tj/ pgs- 
The ohmic time is 



idiff = 2.3 x 10 yrs 



T 



7/2 



i 2 g 2 A e 



(Z) 



Trp'iRl 



and the accretion time is 
facer = 5.1 x 10* yrs 
The ratio of the two is 



5/2 



fdiff =33 Tr_ 

taccr p2 



M- 



(20) 



(21) 



(22) 



/Xp g A ei (Z) 

We have inserted typical values for temperature and density, 
although the exact value of this ratio depends on the T-p 
relation in the atmosphere. For a constant flux atmosphere 
with free-free opacity, T oc p 4 ^ 13 , giving tug /facer oc p~ 3 ^ 13 , 
again weakly dependent on density. 

3.2 Detailed models 

We calculate models of the atmosphere by integrating the 
heat flux equation 

4acT 3 dT 



3k dy 
and the entropy equation 



dF _ cpTm 
dy y 



dlnT 



dlny 



(23) 



(24) 




io- 8 io- 7 io- 6 io- 5 io- 4 io- 3 io- 2 io- 1 

AM (M ) 




Figure 3. (i) Temperature profiles in the white dwarf enve- 
lope for M = 10 — 10 (lower curve) and 10 -9 Mq yr — 1 (up- 
per curve), for a 0.6Mq white dwarf, (ii) Ratio of ohmic time 
to accretion time in the envelope for M = 10 — 10 (lower curve) 
and 10 — 9 Mq yr -1 (upper curve). The composition is solar for 
AM < 1(T 5 Mq, pure He for IO" 5 Mq < AM < IO" 2 Mq, and 
C/O for AM > 1CT 2 M Q . 



where k is the opacity, cp the specific heat at constant pres- 
sure, and the adiabatic index is V a d = (dlnT/dlnP) a d. The 
terms on the right of equation (^) represent the effects of 
compressional heating (Nomoto 1982; Townsley & Bildsten 
2001,2002). 

We integrate equations ( p3| ) and ( |2^ ) from the top of the 
atmosphere to the base, iterating the choice of flux at the 
top until we match the correct flux from the core. We calcu- 
late the opacity, which includes contributions from free-free, 
electron scattering and conduction, as described by Schatz et 
al. (1999). Our calculation of electrical conductivity, which 
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from solar to pure He and from pure He to C/O (eqs. |18| 
and Q give t die /t accl oc 1/ (Z)). At low densities, tdiff/iaccr 
decreases with depth, as expected from equation 

The condition tdm = t accr defines a critical accretion 
rate M c . Inspection of Figure ^ shows that for the accreted 
envelope (AM < IO" 5 Mq), M c « 10~ 10 M Q yr" 1 , and for 
the He layer and the C/O core, an average value is M c « 5 x 
lO" lo A/0 yr" 1 . This agrees well with our analytic estimates 
(compare eqs. |^| and when idiff/tacci = !)• 



4 GLOBAL MODELS OF ACCRETION AND 
DIFFUSION 

In this section, we follow the effects of advection and dif- 
fusion in a simple model of the global magnetic field struc- 
ture. We show that the evolution is very different depending 
on whether M is less than or greater than the critical rate 
M c « (1-5) x 10" 10 M Q yr" 1 found in §3. 



Figure 4. Ohmic diffusion time (soiid fines) and accretion time 
(dotted lines) in the white dwarf envelope for the models of Figure 
| (M = IO" 10 and IO" 9 Mq yr" 1 ). The accretion time (dotted 
lines) is iabelled with the appropriate accretion rate (the accre- 
tion time is smaifer for a higher accretion rate). The diffusion time 
(solid fines) is independent of accretion rate in the degenerate lay- 
ers (the conductivity depends only on density); but is longer for 
higher accretion rates in the non-degenerate layers (the conduc- 
tivity increases with temperature). 



4.1 A simple model of the global field evolution 

We make the following simplifying assumptions: (i) we con- 
sider an axisymmetric, poloidal magnetic field, (ii) we as- 
sume spherical accretion, and (iii) we take the accreted mat- 
ter to have the same composition as the core (equal amounts 
by mass of C and O). Under spherical accretion, each spher- 
ical harmonic component / evolves independently, so that an 
initially dipolar field remains dipolar as accretion proceeds. 
In reality, the accreted matter is channelled onto the polar 
caps in magnetic systems, subsequently spreading over the 
surface of the star. We do not attempt to model this pro- 
cess here, meaning that we must pay careful attention to our 
choice of surface boundary condition. We discuss this issue 
below, but first outline the equation describing the evolution 
of the magnetic field, and how we solve it. 

The induction equation with advection included is 



dB 

dt 



Vx (nxB)-Vx (r/VxB), 



reduces to equation (|3|) for degenerate electrons, and equa- 
tion @ for non-degenerate electrons, is described in the 
Appendix. 

As an illustrative model of the outer layers of an accret- 
ing white dwarf, we take a layer of solar abundance material 
of mass 10" 5 Mq, a pure helium layer of mass 10" 2 Mq, 
and a white dwarf mass of 0.6 Mq, with an equal mixture 
by mass of carbon and oxygen in the core. We choose a lu- 
minosity from the core of 10 -3 Lq. We show the resulting 
temperature profiles for M = 10" 10 and 10" 9 Mq yr" 1 in OT) m terms of the vector potential A, 
Figure 0. For M = 10" 10 Mq yr" 1 , we find a core temper- 
ature 9.4 x 10 6 K, luminosity 4.7 x 10~ 3 Lq and effective 
temperature 13, 600 K. For M = 10~ 9 Mq yr" 1 , we find a 
core temperature 2.0 x 10 7 K, luminosity 6.2 x 10 -2 Lq and 
effective temperature 25, 900 K. These values are in reason- 
able agreement with the detailed models of compressional 
heating recently presented by Townsley & Bildsten (2001, 
2002). A simple estimate of the compressional heating lumi- 
nosity is to write 



(26) 



(27) 



OA 

— = v x ( V x A) - rjV x V x A. 

As previously, we consider an axisymmetric field, so that 
A — A^{r,9)e^. In the absence of diffusion, equation (|27~ 
gives 

/ d 



.dt 



+ vV] (rA^smO) = 



(28) 



L » c P TM « 3.4 x 10" 3 Lq T 7 M. 



(25) 



where we use the heat capacity of an ideal gas cp = 
5kp /2fim p , and take zj « 1 as a mean value in the enve- 
lope. This estimate agrees well with our numerical results. 

The ratio tdiff/taccr is shown in Figure^, and the ohmic 
time and accretion time individually in Figure M. Our numer- 



(WVS; Choudhuri & Konar 2002). The quantity rA^sinfl 
(which labels the magnetic field lines, see §2.2) is advected 
by the poloidal velocity field. 

As before, we separate A$ into radial and angular pieces 
according to equation (0). We take the velocity field to be 
purely radial, v = v r (r,6)e r . Equation ( p7j ) then gives 



ical results agree well with the analytic estimates in §3.1. In 
the degenerate layers (AM > 10 ), tdiff/^accr increases with 
depth (eq. The jumps in tdifi/tacct at AM = 10" 5 M 

and AM = 10 -2 Mq result from the changes in composition 



dRi 
dt 



dRi , 



l(l + l)Ri 



dr 2 



(29) 



This equation governs the evolution of Ri(r,t) under the 
joint action of accretion and diffusion. 
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Figure 5. Time evolution of the magnetic field during accretion of 0.1 Mq onto a 0.6 Mq white dwarf. We assume "vacuum" boundary 
conditions (see discussion in §4.1). The left panel shows M = 10~ 10 Mq yr — 1 , for which the accretion takes 10 9 years; the right panel 
shows M = 10 -9 Mq yr — x , for which the time elapsed is 10 8 years. The curves are spaced by equal increments in central density, which 
increases from s» 3 X 10 6 gem -3 to Ri 6 X 10 e gem -3 . We normalise Ri to have a maximum value Ri = 1 initially, and normalise B r 
and Bg to have a central value of unity initially. 



We evolve equation ( |29|) numerically^. At each time 
step, we include advection by making a new white dwarf 
model (as described in §2.1) with a slightly larger mass, 
but keep Ri fixed for a given fluid element. We then apply 
an implicit Crank-Nicholson scheme for the diffusion term 
(e.g. Press et al. 1992). With accretion switched off, our code 



2 This numerical approach is similar to that of WVS, who in- 
cluded the advection terms in their study of cooling white dwarfs. 
WVS found that the significant contraction which occurs in the 
pre- white dwarf stages of evolution increased the rate of ohmic de- 
cay, by taking the initial n = 1 decay mode and generating higher 
n components. This interesting effect is not significant during the 
accretion process, which involves much smaller changes in white 
dwarf radius. 



accurately follows the predicted decay of the eigenmodes 
found in §2. 

We adopt two different boundary conditions at the sur- 
face. The first is the "vacuum" boundary condition given 
by equation (^). The second boundary condition is a 
"screened" boundary condition Ri = at the surface. In 
the advection step, the boundary condition determines the 
values of Ri assigned to the newly accreted matter. As de- 
scribed above, we take the composition of the accreted ma- 
terial to be equal mass fractions of C and O. At the centre 
of the star, where v r vanishes, we apply equation (p^). 

The "vacuum" boundary condition implies J = out- 
side the star, as in §2, and that each shell of matter is ac- 
creted with the vacuum field at the surface. This boundary 
condition was adopted by Choudhuri & Konar (2002) in 
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Figure 7. Magnetic field lines for the models shown in the left 
panel of Figure 6. We show the initial field structure in the top left 
panel, and the final structure for (clockwise from the top right) 
M = 10~ 9 , 10~ 10 , and lO^ 11 M yr" 1 . The dotted line indicates 
the stellar surface (note the larger radius of the less massive initial 
model). We show field lines with the same value of flux in each 
panel. The dramatic ohmic decay at the lowest accretion rate 
appears as the reduction in the number of field lines. 



their recent study of accreting neutron stars (these authors 
solved eq. |27| with a prescribed 2D velocity field). In our ID 
model, we envisage this boundary condition as approximat- 
ing the case where the accreted material is able to spread 
away from the polar cap without significantly distorting the 
field structure. This is similar to Choudhuri & Konar (2002), 
who put in an outwards radial velocity in the outer layers 
to simulate the effects of buoyancy instabilities. They found 
the evolution of the field was then determined by advection 
in the interior. The "screened" boundary condition Ri — 
assumes screening currents at the surface, similar in spirit to 
the plane parallel models of CZB for accreting neutron stars, 
and to the idea that the magnetic field may be completely 
buried by accretion. 

These two different boundary conditions illustrate the 
range of behaviour that might be expected given a more 
complex flow geometry and better treatment of the outer 
layers, including effects such as buoyancy or interchange in- 
stabilities. However, note that the solution for the interior 
is insensitive to the choice of surface boundary condition. 
In the next section, we present our numerical results. We 
show that the evolution of the field depends upon whether 
the accretion rate is above or below the critical rate M c . 



Figure 8. As Figure [ij, but for the right panel of Figure 6 (screen- 
ing boundary condition). 



4.2 Results 

Figure [| shows the time evolution of a dipole field during 
accretion of 0.1 Mq onto a 0.6Mq white dwarf at two dif- 
ferent accretion rates. We take a "vacuum" boundary con- 
dition, and the initial profile of Ri(r) to be the lowest or- 
der dipole decay mode. We show the time evolution of Ri, 
B r (oc R x /r 2 , eq. §), and B e (oc {l/r){dR 1 /dr), eq. §) 
at equal steps in central density. We normalize Ri so that 
the initial profile has a maximum value Ri = 1, and nor- 
malize B r and Bg so that the initial central value is unity. 
The initial central density is « 3 x 10 6 g cm~ 3 , increasing to 



6 x 10 6 g cm" 3 . The left panel is for M 



10- iu M Q yr" 



in which case the accretion lasts for 10 years, and the right 
panel is for M = 10 -9 Mq yr -1 , in which case the accretion 
lasts for 10 8 years. 

At M = 10~ 10 M Q y^ 1 (left panel of Figure |), the 
magnetic field has time to diffuse into the newly accreted 
material. The overall decrease in Ri is due to ohmic de- 
cay during the 10 9 years of accretion (the lowest order de- 
cay time of the initial model is m = 7.6 billion years). At 
M = 10 -9 Mq yr -1 (right panel of Figure ||), the accreted 
material is added too quickly for significant diffusion, so that 
the surface value of Ri drops as accretion proceeds. The 
central values of B r and Bg increase by a factor of w 1.6. 
This can by understood in terms of the combination of flux 
conservation, which gives r 2 B ^constant, and mass conser- 
vation, which gives r 3 p ^constant. Eliminating r, we find 



B oc 



n 2/3 



giving an increase in B of 1.6 when p increases 



by a factor of 2. (In the thin layer considered by CZB, the 
compression was one-dimensional, giving B oc p in that case; 
see their eq.J13]). 

Figure ^ shows the final magnetic profiles for accre- 
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Figure 6. The profiles of R\, B r , and Bg after the accretion of 0.1 Mq onto a 0.6 Mq white dwarf at different rates. We show results 
for vacuum boundary conditions (left panel) and screening boundary conditions (right panel). The dashed line in each panel shows the 
initial state. The solid lines show accretion rates 10 -8 and 10 — 9 Mq yr _1 , for which mass is accreted too rapidly to become significantly 
magnetised. The dotted lines show accretion rates 10 — 10 and 10 — 11 Mq yr _1 , for which ohmic diffusion has time to magnetise the 
accreted material. Significant ohmic decay of the field occurs at these low rates, because the accretion takes a long time. 



tion at several different rates. We show the profiles for vac- 
uum boundary conditions (left panel) and screened bound- 
ary conditions (right panel). The initial profile in each case 
is shown by the dotted line, and is the lowest order decay 
mode of the initial model (subject to the respective surface 
boundary condition). We show accretion rates M = 10 -11 , 
10" 10 , 1CT 9 and 10~ 8 Mq yr" 1 . At the lowest accretion 
rates, 10 -11 and 10 -10 Mq yr -1 (dashed lines), the magnetic 
field is able to diffuse into the newly accreted layers. The de- 
crease in Ri is due to ohmic decay over the long timescale 
of accretion. For accretion rates 10 -9 and 10 -8 Mq yr -1 
(solid lines), the accreted layers do not have time to become 
significantly magnetized. At M = 10 -8 Mq yr -1 , the final 
profile is almost exactly that given by conserving Ri in each 
fluid element as accretion proceeds. 

In Figures [j] and we show the magnetic field fines 
for the profiles in Figure [| We show field lines with 
the same values of r sin 8 in each panel. The M = 
1CP 11 Mq yr -1 has only two field lines because of the sig- 
nificant ohmic decay that occurs as accretion proceeds. Sim- 
ilarly, at 10~ 9 Mq yr -1 , all of the field lines shown in the 
initial model have been pushed into the interior by the ac- 
cretion flow. 

Figure ^ shows R\, B r , and Bg at the surface as a 
function of mass accreted for a vacuum boundary condi- 



tion. As we discussed §4.1, the quantitative decrease in the 
surface field depends on the chosen surface boundary con- 
dition. However, Figure ^ shows the qualitative effect of ac- 
cretion at different rates. From top to bottom, the accre- 
tion rates are M = KT 10 , KT 11 , and 10" 9 Mq yr" 1 . At 
M = 10 -10 Mq yr -1 , the surface field changes only slightly 
as the newly accreted material is able to become magne- 
tized by diffusion. For M — 10~ n Mq yr" 1 , the surface 
field decreases due to ohmic decay (in this case the accre- 
tion takes 10 10 years, comparable to the decay time). For 
M = 10~ 9 Mq yr -1 , the change is directly due to accretion. 



5 IMPLICATIONS FOR OBSERVED SYSTEMS 

In the previous sections, we have shown that, under the as- 
sumption that the white dwarf retains the mass it accretes, 
the surface magnetic field may be significantly reduced for 
accretion rates greater than the critical value M c m (1- 
5) x KT 10 M @ yr' 1 . We now conclude by discussing the 
implications of our results for observed systems. We show 
that there are several interesting consequences of abandon- 
ing the assumption that the magnetic fields of accreting 
white dwarfs do not change with time. 
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Figure 9. For the vacuum boundary condition solutions, we 
show surface values of R\ (dotted line), B r (solid line), and 
Bg (dashed line) relative to their initial values as a function 
white dwarf mass for (top to bottom) M = 10 — 10 , 10 — , and 
10" 9 Mq yr" 1 . 
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Figure 10. The orbital periods of AM Her systems and inter- 
mediate polars. We show orbital periods for 53 AM Her systems 
from Wickramasinghe & Ferrario (2000), and for 20 intermediate 
polars from the catalogue of Ritter & Kolb (1998). We exclude 
the intermediate polar GK Per which has P or b = 2 days and an 
evolved companion. 



5.1 An evolutionary connection between 
AM Hers and intermediate polars? 

Accreting white dwarfs fall into two classes. Those in AM 
Her systems rotate synchronously with the binary orbital 
period, and have fields 10 7 -3 x 10 8 G (WF). Intermediate 
polars (see review by Patterson 1994), in which the white 
dwarf rotates asynchronously, are believed to contain white 
dwarfs with fields ~ 10 5 -10 7 G, although most have not 
been measured directly (the lower limit here is the field 
strength needed to disrupt the accretion flow before it hits 
the white dwarf surface). The asynchronous rotation of the 
white dwarfs in IPs lead to early suggestions that the white 
dwarf magnetic fields are lower than in AM Hers (for ex- 
ample, Lamb & Patterson 1983 estimated B ~ 10 6 G). In 
addition, whereas strong optical polarization is seen in AM 
Hers, little or no polarization is seen in IPs, also suggest- 
ing weaker magnetic fields. Recently, the lack of Zeeman 
splitting in spectroscopic observations of the white dwarf in 
V709 Cas (Bonnet-Bidaud et al. 2001) has constrained the 
magnetic field to be < 10 7 G. 

The two classes show an interesting difference in ac- 
cretion rates. AM Her systems accrete at low rates, M as 
5x 10" 11 M yr _1 (Chanmugam, Ray, & Singh 1991; Warner 
1995), whereas the IPs accrete rapidly, Warner (1995) esti- 
mates M « (0.2-4) x 10" 9 M e yr" 1 from X-ray fluxes of 
ten systems. Wickramasinghe & Wu (1994) propose that 
the low accretion rate of AM Hers is a direct consequence of 
the strong magnetic field of the white dwarf, which inhibits 
magnetic braking by reducing the number of open field lines 
along which the magnetized wind from the secondary may 
flow (for example, see Li, Wickramasinghe, & Wu 1995). 

Figure [To] shows the orbital period distributions of 
known AM Hers and IPs. Most AM Her systems have 
forb ^ 3 hours, whereas most IPs have P or b « 3-6 hours. 
This difference has led many to suggest an evolutionary link 
between the two classes, namely that the asynchronous IPs 
are the progenitors of the synchronous AM Her systems. 
For a magnetic field B ~ 10 7 G, the magnetospheric ra- 
dius equals the orbital separation at P or b ~ 4 hrs, suggest- 
ing a natural dividing line between synchronous and asyn- 
chronous systems (Chanmugam & Ray 1984; King, Frank, 
& Ritter 1985; Hameury et al. 1987). However, an evolution- 
ary connection is difficult if the white dwarf magnetic field 
is constant with time (e.g., Wu & Wickramasinghe 1993). It 
requires either (i) the IPs and AM Hers have similar mag- 
netic field strengths, but the optical polarization is somehow 
suppressed in asynchronous rotators, or (ii) the IPs have 
weaker fields than the AM Her systems, implying a yet to 
be detected population of asynchronous AM Her progenitors 
(King & Lasota 1991). 

It is striking that the mean accretion rate of the IPs 
is > M c , whereas the AM Her systems accrete at rates < 
M c . Thus, given our results, a natural suggestion is that the 
magnetic fields of the intermediate polars are low because 
of the direct effects of the accretion flow. If so, it is possible 
for an intermediate polar to evolve into an AM Her system 
if the mass transfer rate to drops from > M c to < M c at 
an orbital period w (3-4) hours (between the observed AM 
Her and intermediate polar populations). 

The reason for such a change in accretion rate is not 
clear, but is reminiscent of the need for a rapid drop in 
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mass transfer rate to explain the period gap (P or b ~ 2-3 
hours) in non-magnetic CVs. Whether or not the magnetic 
CVs exhibit a period gap is unclear, there are certainly a 
higher fraction of magnetic systems in the gap than for non- 
magnetic systems (King 1994; Kolb 1995; Wheatley 1995). It 
is however possible that a similar mechanism is operating in 
both cases, perhaps disrupted magnetic braking, as preferred 
for non-magnetic systems (Rappaport, Verbunt & Joss 1983; 
Howell, Nelson, & Rappaport 2001). 

There is an alternative scenario along the same lines as 
that of Wickramasinghe & Wu (1994). The ROSAT detec- 
tion of IPs with soft X-ray spectra similar to weaker field AM 
Hers (Haberl & Motch 1995) has led to suggestions of two 
populations of IPs. A possibility is that those IPs with fields 
close to 10 7 G come into synchronous rotation at P or b ~ 4 
hours, which causes a drop in M as suggested by Wickra- 
masinghe & Wu (1994), allowing the field to subsequently 
grow to values ~ 10 s G typical of AM Hers. A problem with 
this picture is what happens to the weaker field IPs (i.e., 
those with traditional hard X-ray spectra), since very few of 
them are seen at short orbital periods. 

The timescale for reemergence of the field is set by 
the ohmic diffusion time at the base of the accreted layer. 
Figure 4 and equation (16) show that this is w 3 x 
10 8 years (AM/0. 1M©) 7/5 where AM is the accreted mass. 
This is similar to the time for a non-magnetic CV to cross 
the period gap (Howell et al. 2001). Thus if the drop in mass 
accretion rate causes the secondary to lose contact with the 
Roche lobe, there may be time for the field to reemerge be- 
fore contact is resumed. The impact of any gap on the ob- 
served orbital period distribution depends on how strongly 
the transition orbital period depends on mass transfer rate 
and magnetic field. With or without a period gap, the lack 
of known synchronous low field systems could constrain the 
timescale for reemergence and therefore accreted massj^]. 

A detailed test of either scenario requires comparison of 
evolutionary studies and the observed populations of IPs and 
AM Hers. We note that whereas there is a clear difference 
in orbital periods and magnetic field strengths between AM 
Hers and IPs as classes of objects, no correlation is observed 
between magnetic field strength and orbital period amongst 
AM Her systems only (WF). 

5.2 Where are the 10 9 G accreting WDs? 

The highest measured magnetic field in an accreting sys- 
tem is 2.3 x 10 8 G for AR Uma, whereas the highest ob- 
served fields in isolated WDs are close to 10 9 G. Whether 
or not this difference is significant remains to be seen. For 
example, WF argue that selection effects hide the very high 
field accreting WDs. There have been suggestions as to why 
the observed upper limit might differ. Wickramasinghe & 
Wu (1994) argue that because the high field systems come 
into synchronization at wider separation, magnetic braking 
is suppressed earlier, and they do not have time to evolve 
via gravitational radiation into contact (thus they predict a 
population of high field, detached binaries) . 

We would like to point out that some difference would 
be expected simply because of the different ohmic decay 

3 We thank the referee, J.-M. Hameury, for bringing up this issue. 



properties of the two populations. The liquid accreting white 
dwarfs decay at a faster rate than the isolated white dwarfs, 
which form a solid core after ~ 10 9 years of cooling, effec- 
tively switching off ohmic decay. To find the maximum dif- 
ference between the two populations, we assume the isolated 
white dwarfs undergo no magnetic field decay. The contrast 
in field strengths is then exp(t/V), where r is the decay time 
for the accreting population. For r = 10 10 years, this gives 
a factor of 2 after t = 7 billion years, or a factor of 3 after 
t — 11 billion years. This effect alone may not explain the 
observed factor of w 3, but should be taken into account 
when comparing the maximum field strengths of the two 
populations. 

5.3 Type la progenitor systems 

Currently favoured models for Type la supernovae progen- 
itors are rapidly accreting (M>10 -7 Mq yr _1 ) massive 
white dwarfs, which are able to burn the accreted hydrogen 
and helium steadily and grow to near-Chandrasekhar mass. 
Examples are the white dwarfs in supersoft X-ray sources 
and symbiotic binaries. If accretion is acting to suppress the 
surface magnetic field in these sources as we propose for 
intermediate polars, we would not expect to find strongly 
magnetic white dwarfs in these systems. There is little in- 
formation about the incidence of magnetism in these rapidly 
accreting systems. In a survey of 35 symbiotics, Sokoloski & 
Bildsten (1999) (see also Sokoloski, Bildsten, & Ho 2001) 
discovered a magnetic WD in the symbiotic binary Z And 
rotating with a 28 minute spin period. If the white dwarf 
is in spin equilibrium with an accretion disk at the mag- 
netospheric radius, the implied magnetic field strength is 
6 x 10 6 G, typical of an intermediate polar. XMM obser- 
vations of M31 revealed a supersoft X-ray source with 865s 
pulsations, interpreted as the spin of a magnetized white 
dwarf (Osborne et al. 2001). A similar estimate for the mag- 
netic field gives « 10 7 G (King, Osborne, & Schenker 2002). 
It remains to be seen how common systems like these are 
amongst rapidly accreting sources. 

As the central density increases due to accretion, the in- 
ternal field is amplified by compression (B oc p 2 ^ 3 ; see §4.2). 
This amplification is shown in Figure 5 for accretion onto a 
0.6 Mq white dwarf, and would be even greater for a mas- 
sive white dwarf. The density increases especially quickly as 
you approach Chandrasekhar mass and the electrons become 
more relativistic, as pointed out by work on compressional 
heating (Nomoto 1982; Townsley & Bildsten 2001, 2002). 
An important problem is to understand the role played by 
such an amplified magnetic field in the approach to and evo- 
lution after ignition of a Type la supernova. For example, 
Ghezzi et al. (2001), recently pointed out that a large scale 
magnetic field of 10 8 -10 9 G can generate asymmetries in the 
thermonuclear flame front that engulfs the white dwarf dur- 
ing the supernova. Our results show that the prior evolution 
of the field during accretion should be taken into account. 

5.4 Complexity as a result of accretion 

Both isolated and accreting white dwarfs often show a 
complex magnetic field geometry, requiring a significant 
quadrupole component or an offset dipole (WF). In §2, we 
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found decay times of (4-6) x 10 9 years for the I = 2 decay 
mode of liquid white dwarfs. Thus one possibility is that 
the quadrupole field is fossil, requiring that the dipole and 
quadrupole fields had comparable magnitudes initially. 

In isolated systems, the liquid decay time is relevant 
only for the first ~ 10 9 yrs, after which the formation of a 
solid core increases the ohmic time. Muslimov et al. (1995) 
studied the Hall effect in isolated white dwarfs as a possi- 
ble mechanism for generating field complexity as the white 
dwarf cools. However, the Hall effect is not expected to play 
a significant role for a liquid interior (Muslimov et al. 1995). 

Another possibility is that the accretion flow directly 
leads to complexity of the field. The spreading of material 
away from the polar cap (which we have not included in our 
models in this paper) could induce higher order components 
of the surface field. An interesting observation is that the 
magnetic pole undergoing most accretion in AM Her systems 
is the pole with the weaker magnetic field in all cases with 
magnetic field measurements for both poles (WF). Wickra- 
masinghe & Wu (1991) suggest that this is a result of the 
role of the quadrupole component during synchronisation. 
An alternative suggested by our results is that after syn- 
chronisation occurs, the pole pointing towards the compan- 
ion undergoes more rapid accretion which reduces the field 
strength at that pole, giving rise to the observed asymmetry. 
Even though the global rate for AM Hers is < M c , the local 
accretion rate in the layers for which the accreted material 
is still confined to the polar cap will be greater than the crit- 
ical rate. Detailed models of the spreading of material away 
from the polar cap are required to investigate this further. 



6 CONCLUSIONS 

We have studied the evolution of the magnetic field in an 
accreting white dwarf. Previous studies of ohmic decay in 
isolated white dwarfs showed that the ohmic decay time 
is always longer than the cooling time because the white 
dwarf develops a solid core. In §2, we calculated the ohmic 
decay times for the lowest order modes of accreting white 
dwarfs, which have a liquid interior because of compressional 
heating by accretion. We found that the lowest order ohmic 
decay time is (8-12) x 10 9 yrs for a dipole field, and (4- 
6) x 10 9 yrs for a quadrupole field (see Table 1), in good 
agreement with the earlier calculations of WVS. The differ- 
ence in ohmic decay times between isolated and accreting 
white dwarfs should be taken into account when compar- 
ing the maximum field strengths of the two populations. 
In addition, the decay timescale for the quadrupole is long 
enough that observed quadrupole components (see WF for a 
summary of the observations) in both accreting and isolated 
white dwarfs may be fossil if the quadrupole is initially of 
similar strength to the dipole component. 

In §3, we compared the timescale for ohmic decay with 
accretion, and showed that accretion occurs more rapidly 
than ohmic diffusion for accretion rates greater than the 
critical rate M c « (1-5) x 1CT 10 M yr -1 . In §4, we calcu- 
lated the time evolution of the magnetic field as a function 
of accretion rate, assuming the white dwarf mass increases 
with time. For a simplified field and accretion geometry (an 
axisymmetric poloidal magnetic field and spherical accre- 
tion), we found that accretion at rates M > M c leads to a 



reduction in the field strength at the surface, as the field is 
advected into the interior by the accretion flow. 

The main conclusion of this paper is that, due to the 
direct action of accretion, significant changes in the surface 
magnetic field of an accreting white dwarf could occur dur- 
ing its accretion lifetime. In §5, we showed that accretion 
induced magnetic field evolution may explain several fea- 
tures of observed systems. Most striking is that the strongly 
magnetic (B ~ 10 7 -3 x 10 8 G) AM Her syste ms have a 
mean accretion rate w 5 x 1CP 11 M© yr -1 < M c whereas 
the weakly magnetic (B ~ 10 5 -10 7 G) IPs accrete at rates 
~ 10~ 9 M e yr _1 > M c . This raises the possibility that the 
white dwarfs in IPs have subsurface fields as strong as those 
in AM Hers. If so, this allows for an evolutionary connec- 
tion between the long orbital period (P or b ~ 3-6 hours) IPs 
and short orbital period AM Hers (P or b < 3 hours), requiring 
that evolution in orbital period somehow causes a drop in ac- 
cretion rate (for example, by disrupted magnetic braking as 
postulated for non-magnetic CVs; see §5.1 for further discus- 
sion). The magnetic field would then reemerge on a timescale 
set by ohmic diffusion, « 3 x 10 s yrs (AM/0.1 M ) 7/5 , 
where AM is the accreted mass. Evolutionary calculations 
are necessary to compare this scenario with the properties 
of observed systems. We did not consider non-magnetic ac- 
creting systems in §5; however, it is possible that many 
non-magnetic white dwarfs (meaning B < 10 5 G so that the 
accretion flow is not disrupted) have submerged magnetic 
fields if they are accreting at rates > M c . 

The models presented in this paper are simplified, and 
many outstanding theoretical questions remain. The ma- 
jor uncertainty in the evolution of accreting white dwarfs 
is whether the white dwarf mass increases or decreases with 
time. Even if the mass is increasing with time as we have as- 
sumed here, classical nova explosions eject some mass from 
the system; the effective mass accretion rate is therefore less 
than the mass transfer rate from the secondary by some fac- 
tor, introducing an extra uncertainty. If classical novae eject 
more mass than accreted prior to the thermal runaway, ex- 
cavation of the white dwarf core leads to a decreasing mass 
with time. This raises the possibility that the surface mag- 
netic field could grow with time if underlying magnetic field 
is exposed as the mass decreases. This was originally sug- 
gested as a source for the magnetic field of DQ Her by Lamb 
(1974). Interestingly, this alternative is also consistent with 
an evolutionary link between AM Hers and intermediate po- 
lars, since the magnetic field would increase as the orbital 
period decreases and successive nova explosions occur. The 
evolution of the magnetic field under the action of mass loss 
is left to a future study. The result is presumably somewhat 
dependent on the initial choice for the current distribution, 
and differences in the thermal evolution will have to be taken 
into account. 

It is important to understand the spreading of matter 
away from the polar cap. This topic has received little at- 
tention, and is not well-understood. Studies of the polar cap 
of accreting neutron stars (Hameury et al. 1983; Brown & 
Bildsten 1998; Litwin, Brown, & Rosner 2001) and white 
dwarfs (Livio 1983; Hameury & Lasota 1985) suggest that 
spreading occurs once the sideways force due to the hydro- 
static overpressure overcomes the confining magnetic ten- 
sion, or sooner if instabilites play a role (Litwin et al. 2001). 
Understanding the spreading process, and its affect on the 
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magnetic field geometry in the outermost layers of the star 
is directly relevant for interpreting observations of polar cap 
fields of AM Hers. A related issue which remains to be ex- 
plored is the stability of the magnetic field configurations 
found in §4. For example, stability considerations demand a 
toroidal field component which is not included in our simple 
models (see discussion in Mestel 1999). 

Finally, we have not discussed the effect of white dwarf 
rotation. As mentioned in the introduction, King (1985), ap- 
plying the results of Moss (1979) to accreting white dwarfs, 
suggested that meridional currents could be responsible for 
the submergence of magnetic flux in the outer layers of ro- 
tating magnetic white dwarfs. The aim was to explain the 
narrow range of field strengths in AM Hers. The effect of 
meridional flows in accreting white dwarfs deserves further 
investigation. 
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APPENDIX A: ELECTRICAL CONDUCTIVITY 
FOR ARBITRARY DEGENERACY 

In the white dwarf envelope, the electrons go from being non- 
degenerate to degenerate. Hubbard & Lampe (1969) cal- 
culated the thermal conductivity for arbitrary degeneracy, 
presenting their results in tabulated form. Unfortunately, 
other calculations of conductivity which present their re- 
sults in terms of fitting formulae, such as those of Yakovlev 
& Urpin (1980) and Itoh et al. (1983) are for degenerate 
electrons only (for thermal conductivity, this is the regime 
of interest, since radiative heat transport dominates at low 
densities). In their study of ohmic decay in cooling white 
dwarfs, WVS adopted an interpolation scheme between cal- 
culations in the degenerate and non-degenerate regimes. We 
have taken a different approach, which we describe below. 
After our calculations were completed, we learned of the new 
conductivity calculations by Potekhin (1999) and Potekhin 
et al. (1999) which treat the regime of intermediate degen- 
eracy] 4 ]. Comparing with these calculations, we find that our 
method agrees to within 20% for degenerate electrons and to 



within a factor of 2 for non-degenerate electrons, adequate 
for our purposes. 

For non-relativistic electrons, the density and electrical 
conductivity may be written 

_ V2(m e k B T) 3/2 



-Fi/2 (a) 



and 



2m e (k B T) 3 
a — — F2 (a) . 



(Al) 



(A2) 



where a = —{Ep — m e c 2 ) /k B T ', and the Fermi integrals 
F n (rj) are given by 



F n (a) 



x n dx 



1 + exp(x + a) ' 



(A3) 



In equation (A2), we assume the Coulomb logarithm A e 



varies slowly enough with electron energy that it may be 
taken out of the integral. For degenerate electrons, a <C 
— 1, equation (A2) reduces to the small x (non-relativistic) 
limit of e quat ion ([|). For non-degenerate electrons, a> 1, 
equation (A2) reduces to equation ( |l9| ) (Spitzer 1962). 

Our approach in this paper is to evaluate a using equa- 
tion ( |A2|) , with Aei a function of density, and with a correc- 
tion for relativistic effects. We evaluate the Fermi integral 
F2(a) by direct numerical integration. The electron nu mber 
density is used to determine a, by inverting equation ( Al) 
using the fit of Antia (1993) to the inverse of F%/2- 

We calculate the Coulomb logarithm as a function of 
density by interpolating between the non-degenerate result 
of Spitzer (1962) and the degenerate results of Yakovlev 
& Urpin (1980). We adopt the expression of Yakovlev & 
Urpin (1980) for the Coulomb logarithm, but rewrite it 
in terms of the electron momentum, choosing x in such a 
way that in the non-degenerate limit, we obtain Spitzer's 
formula. The Coulomb logarithm is A e i = A^ — v 2 /2c 2 , 
where A™ = ln(r max /r min ) (YU). Here r min and r max are 
the limits of the integral over impact parameters. For tem- 
peratures greater than 4.2 x 10 K (Spitzer 1962), r m i n is 
set by the de Broglie wavelength of the electrons, we take 
r m i n ~ h/2p e . The cutoff r max is set either by the Debye 
length, given by r 2 D = m p k B T / Aite 2 p ^ YiZ 2 , or the inte- 
rion spacing a = (3/47rni) 1//3 , whichever is larger. We obtain 



= In 
= In 



/ 2 X 1 / 3 m e c x /3 3\ 

\9^) ~(j2n t y/ a \ f + 2) 

_/3 3X 1 / 2 

1/3 Vr + 2) 



3X 1 / 2 



127 x 



where 



p 1/3 (E^) 
Enaf _ Qll fE i Y i zf\ (p 5 E^) 1/3 



k B Ta j 



T 8 



(A4) 



(A5) 



(Hubbard & Lampe 1969; F = (l/3)(a/r D ) 2 ), and where 
x = (x{ + xj) 1/2 , with x 2 = (3k B T/m e c 2 ) 1/2 = 0.22 T 8 1/2 . 

Finally, we include a cor rection for relativistic effects. 
Because we use equation ( |Al[ ) to solve for a, we obtain the 



4 We thank Alexander Potekhin for bringing these results 
to our attention. Details of these calculations and com- 



nntnr nodes to calculate conductivities can he found at 



tittp : //www. ioffe.rssi . ru/astro/ conduct /index . html 
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non-r elat ivistic value for a in the degenerate limit, and equa- 
tion (A2) reduces to the small x limit of equation (^). To 
include relativistic effects, we divide equation ( A2) by a fac- 



tor 1 + x , so that we recover equation (g) in full in the 
degenerate limit. 

For the T = models of §2 and §4, we write out A^ 
as a sum of logarithmic terms, but drop the term containing 
the temperature T, which contributes < 10% to the final 
value of A ei . 
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